Orbital ordering in LaMnO^: Electron-lattice versus electron-electron interactions 
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The relative importance of electron-lattice (e-1) and electron-electron (e-e) interactions in ordering 
orbitals in LaMnOa is systematically examined within the LDA+t/ approximation of density func- 
tional theory. A realistic effective Hamiltonian is derived from novel Wannier state analysis of the 
electronic structure. Surprisingly, e-1 interaction (~ 0.9 eV) alone is found insufficient to stabilize 
the orbital ordered state. On the other hand, e-e interaction (~ 1.7 eV) not only induces orbital 
ordering, but also greatly facilitates the Jahn- Teller distortion via enhanced localization. Further 
experimental means to quantify the competition between these two mechanisms are proposed. 
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Study of perovskite manganites has been one of the 
main focuses of recent research in condensed matter 
physics, not only because of their great potentials in tech- 
nological applications related to colossal magnetoresis- 
tance (CMR) in Lai-xCa^MnOa, but also because these 
strongly correlated electron materials (SCEM) are ideally 
instrumental to the understanding of the complex inter- 
play of the charge, spin, orbital, and lattice degrees of 
freedom that leads to abundant fascinating phenomena, 
including long range order in all the above channels 0. 

The unusual orbital degree of freedom in the mangan- 
ites, which is the focus of our study, originates from the 
singly occupied degenerate e g states (d z 2 and d x 2_ y 2) 
of the Mn 3+ 3d electrons in the high-spin configuration 
(^■2g e g) due to the ligand-ficld splitting and strong Hund's 
coupling. This orbital degeneracy makes the Mn 3+ ion 
Jahn- Teller (JT) active: the degeneracy can be split via 
biaxial distortion of the surrounding oxygen octahedron. 

Currently, one of the critical questions on the man- 
ganites is the interplay of electron-lattice (e-1) and e g 
electron-electron (e-e) interactions 0, 0, 0, S S 13 • 
Indeed, carrier mobility 0] and magnetism Ul^, both 
essential to the unresolved mechanism of CMR in the 
manganites 0, 0, 0, 0] , respond to these two inter- 
actions in a substantially different manner (despite some 
other aspects reacting similarly). However, even for 
the simplest parent compound, LaMnOa (which presents 
prototype orbital order (00) and strongest JT-distortion 
in the family) . various spectral measurements to date 
[3 QH El Hi GJ] have left two possible mechanisms 
of the OO in dispute: the cooperative JT e-1 effect EH 
and the e-e superexchange effect (2tJ . To our knowledge, 
there is no clear experimental evidence to support one 
over the other, due to lack of "signatures" in distinguish- 
ing these two mechanisms. This leads to great confu- 
sions in the field: for example, whether a new type of 
elementary excitations called orbiton has been observed 
in Raman scattering spectroscopy |l7tll8l |. and more gen- 
erally which interaction dominates localization of the e g 
electrons and thus facilitates the CMR effects upon dop- 



ing LaMn0 3 

It is thus important and timely to 
discern the real roles of e-1 and e-e interactions. 

In this Letter, the electronic structure of the prototype 
LaMnOa is systematically analyzed, aiming to quantify 
the relative importance of e-e and e-1 interactions in or- 
dering the orbitals. A realistic effective Hamiltonian for 
the low energy e g states is derived from a novel Wan- 
nier function-based scheme for general SCEM, with the 
key effective e-e interaction U e s ~ 1.7 eV, JT splitting 
Ajt — 0.9 eV, and octahedral-tilting induced tetrago- 
nal crystal field E z ~ 0.12 eV. Surprisingly, c-1 interac- 
tion alone is found insufficient to stabilize OO. On the 
other hand, e-e interaction not only induces OO, but also 
greatly facilitates the JT distortion by strongly localiz- 
ing the electrons. The present results provide new insight 
into OO in LaMnOs, and place stringent constraints on 
any realistic theories of excitations and CMR in the man- 
ganites. Furthermore, our analysis indicates certain com- 
petition between mechanisms, allowing direct experimen- 
tal determination of their relative strengths. 

Aiming at a careful determination of the relevant 
mechanisms without assuming a particular dominant in- 
teraction, the following novel three-step scheme designed 
for general SCEM is employed, with each step leading to 
additional insights into essential interactions and more 
accurate evaluation of their strengths: (i) A systematic 
analysis of the electronic structure within the LDA+£7 
method [12, against various constraints (e.g.: lattice 
distortion), giving a rough estimate of the relevant e-1 
and e-e effects; (ii) Construction of the low-energy Wan- 
nier states 0, |21j , leading to signatures of competition 
between e-1 and e-e interactions, additional tetragonal 
field effects, and a local representation of the LDA+C7 
Hamiltonian, ijLDA+^j (iii) Determination of the effec- 
tive many-body Hamiltonian via a self-consistent map- 
ping to 7^ LDA +c / ; providing a quantitative evaluation of 
the interactions and deep microscopic insights. 

(i) Analysis of the LDA + U results. The results of our 
systematic study of the electronic structure of LaMnOs 
are summarized in Fig. ^ where the total energy gain per 
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FIG. 1: Total energy gain per formula unit via JT distortion of 
the MnC>6 octahedra, expressed in Q, = (Qf ,Qf) = (\/2(2 — 
s), ^2/3(2771— l—s)), with I, m, and s being the long, medium, 
and short Mn-0 bond lengths, respectively. For cooperative 
JT distortions in LaMn0 3 , Qi = (Q cos 6q , e iq ' R< Q sin 6q ) at 
the i-th Mn site located at R;; in low-temperature phases, 
Q = Q* = 0.4 A, 9 q = 8q = 109°, and q= (ir,7r,0) ~ 



formula unit (f.u.), Ejt, as a function of cooperative JT 
distortion vectors Qi (defined in the caption) is shown 
for a wide range of U [23L [24| . Notice that for realistic 
U = 8 cV, the crystal structure is stabilized at Q* = 
0.4 A [Fig. [Ha)] and 6* Q = 109° [Fig. □»], in excellent 
agreement with experiments [27J, supporting the good 
quality of the LDA+J7 approximation for this system. 

Surprisingly, e-1 interaction alone is found insufficient 
to stabilize the orbital ordered insulating phase: for 
U = eV, the system stabilizes in a metallic state at 
Q = 0.8Q* with small Sjt = —27 meV, despite being 
weakly insulating at Q — Q* with a small gap of 0.1 
eV |2jj. On the other hand, Ejt dramatically increases 
in magnitude to —215 meV as U increases to 8 eV, in- 
dicating that the electron localization induced by the e-e 
interaction greatly facilitates the JT instability — in fact, 
OO can be stabilized even without JT distortion [22^ . 
Indeed, as shown in Fig. ^a), fitting the data points to 
the JT picture {E JT ~ -\gQ + \KQ 2 [l]|) is excellent 
for U > 4 eV but unsatisfactory for small U, suggesting 
that only with U > 4 eV the e g electrons are well local- 
ized as assumed in the JT picture. For realistic U = 8 
eV, the JT splitting (~ 4E JT at Q*) is thus A JT ~ 0.9 
cV, comparable to ~ 0.8 eV extrapolated from spectral 
cllipsometry [lillltij. 

It is extremely important to distinguish carefully the 
JT distortion from other lattice distortions, namely octa- 
hedral tilting and its associated octahedral distortion of 
GdFe03-type. The latter is found negligible in LaMnOs, 
as our study of LaFe03 (a JT-inactive counterpart of 
LaMn03 |2Jj) produces negligible octahedral distortion 
with total energy gain of merely —3 meV. In contrast, 
octahedral tilting is found to give large energy gain of 
—403 meV/f.u. Clearly, a correct assessment of Ejt 
should exclude such a JT-distortion-unrclated contribu- 
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FIG. 2: Left panel: LDA+C7 (8 eV) band structure of 
LaMnOs in the real crystal structure, with circles reflect- 
ing the weights of the Mn d z n (black) and d x 2_ y 2 (blue) 
symmetries. M=(tv, 0, 0), r=(0, 0, 0), R=(ir/2, tt/2, it/2), and 
A=(tt,0,tt/2). Solid green and red lines show the dispersions 
of OO-relevant spin-majority e g Wannier states, illustrated in 
right panels (see text). 



tion. Therefore, our systematic study was performed 
with the experimental octahedral tilting angle (0tiit=16°, 
nearly unchanged upon JT distortion |26l l27| ). 

In addition to localizing e g states, e-e interaction 
plays other crucial roles, as clearly demonstrated in the 
LBA+U (8 eV) band structure (Fig.0 left panel): The 
~ 2.6 cV splitting between the spin-majority occupied 
and unoccupied e g bands near the Fermi energy is too 
large to be accounted for with the estimated JT split- 
ting (~ 0.9 eV), indicating an effective on-site repulsion 
U c g ~ 1.7 eV. Note that the U e s relevant to OO is to be 
distinguished from the "bare" U acting between atomic 
Mn 3d states, as Z7 e ff acts only between the low-energy e g 
Wannier states (WSs, discussed below), and thus includes 
effects of additional screening and slight derealization 
via hybridization that weaken the "bare" repulsion. Also 
note that the considerable amount of e s -character near 
the bottom of the oxygen 2p bands ([—8, —6] eV) is not 
very relevant to OO, as such feature has its origin in 
strong hybridization with the oxygen 2p orbitals, inde- 
pendent of ordering of the orbitals. Though it does imply 
the charge-transfer nature of the manganites in general. 

(ii) Construction of WSs. To proceed with more quan- 
titative evaluation of the above effects and to identify 
other relevant mechanisms, a well-defined local (Wan- 
nier) representation of the low- energy e g states is neces- 
sary for further theoretical formulation. To this end, our 
previously developed energy-resolved symmetry-specific 
WS construction [lj,|2lj is extended to allow mixed sym- 
metry with a constraint search for maximal localization 
[29| . The resulting orbital ordered occupied ([-2.5,0] eV) 
and unoccupied WSs ([0,2.5] eV) are illustrated in Fig. [21 
(middle panel), from which the staggered ordering of the 
orbitals is apparent, as well as the considerable weight 
at the oxygen sites due to strong p-d hybridization. Also 
given in Fig. [21 are band dispersions (solid lines) corre- 
sponding to these WSs, obtained by diagonalizing the 
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FIG. 3: (a) 0t as a function of JT distortion, (b) Compari- 
son of results from Eq. 10 (lines) and LDA+C/: 
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LDA+C7 Hamiltonian in Wanner representation, 



H 



LDA+C/ 
' jm' ,im 



(jm'\H 



LDA+C/i 



(1) 



where \im) denote the m-th WS at site i. For convenient 
theoretical formulation, a second set of WSs with pure 
d 3z 2_ r 2 and d x 2_ y 2 symmetries are built from unitary 
transformation ('remixing', c.f. Fig.JSJ of the first set of 
WSs. The resulting "conventional" WSs form a realistic 
basis for our theoretical modeling below. 

An important observation emerges from the occupied 
WSs in Fig. [21 They do not possess the d^ x 2_ r 2 and 
d 3y 2_ r 2 symmetry expected from the JT picture |l9j . 
Indeed, the occupied states \occ.) = cos^-\3z 2 — r 2 ) ± 
sin^lx 2 - w 2 ) do not follow 6 T = Oq = 120° from the 
JT model |l9|. As it will become more apparent with 
the realistic Hamiltonian below, while e-e and e-1 inter- 
actions induce cooperatively the staggered OO in this 
system, different mechanisms actually compete in deter- 
mining the 'orbital mixing angle', 6t- Specifically, the 
purely electronic (superexchange) mechanism |20j favors 
9t = 90° for the cubic perovskite structure 0, Jj, which 
was confirmed by our LDA+C7 calculations (not shown). 
This competition makes 9^ a sensitive measure of the 
relative importance of leading mechanisms for OO. 

An additional relevant mechanism can now be clearly 
observed from Fig.|3Ja), in which 9?, evaluated from the 
above transformation between the above two sets of WSs, 
is plotted against various magnitude of JT distortion Q 
with experimental 9q. In the absence of the JT effect 
(Q = 0), 6t = 76° is much smaller than 90° expected 
from the superexchange effect 0, 0|- This reflects the 
importance of the tetragonal crystal field, E z , yielded 
mainly from the octahedral-titling in the real structure 
0, as E z favors 6t = 0° instead. Notice also from 
Fig. Efa) that 9t grows as Q increases, reflecting a size- 
able e-1 interaction. Nevertheless, even at optimal (ex- 
perimental) Q, 9t is still way below 9q, confirming that 
the JT effect is far from being dominant. 



(iii) Mapping the Effective Hamiltonian. Based on the 
above systematic analysis, the concrete physics of the 
low-energy spin-majority e g states can be described by 
an effective Hamiltonian including U e s, Ajt, and E z : 

RCS = J- t%'dl y d n + U cS Y,n tl n,l 

- g Ti • Qi - E z £ Tf + K(Q), (2) 

i i 

where 7 and 7' refer to the conventional WSs, denoting 
I T) = |3z 2 -r 2 ) and | |) = \y 2 -x 2 ). T, = (I?,2f) with 
- (4 T d, T - dld n )/2 and T? = (d^ + 4^/2 
are pseudo-spin operators. The in-plane nearest neigh- 



tjj = t'lj = ^fy^t/A (the =p sign distinguishes hopping 
along the x and y directions), while out-of-plane hop- 
pings are strongly suppressed in the A-type antifcrro- 
magnetic ground state due to the double-exchange effect 
Q 9 is the JT coupling constant and K{Q) = 
1 12i Qj 1 Q» + ■ ■ ■ consists of lattice energies. 

To properly map out these parameters from the 
LDA+f/ results, an unambiguous approach is developed 
in this study. Employing the fact that strong local e-e 
interaction is approximated in LDA+J7 in an effective 
HF manner [22|, a proper connection can be made on 
the same WS basis by matching H^5 ^ U with the self- 
consistent HF expression of H ctf : 
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bor hopping integrals are t\- = 3t/4, t\- = t/4, and 
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t 11 ' d} ,d- 



Yj • Bj + Eq . 



(3) 



Here the effective field = (Bf,Bf), where Bf = 
2U cS (T,?)+E z +gQ! and Bf = 2U cft (Tf ) +gQ* . {T?) - 
T z and (Tf) = e^ Ri T x with q = (7r,7r,0) are the com- 
ponents of the "pseudo-spin density wave" mean field 
(T;), which must be determined self-consistently and is 
not necessarily parallel to Q;. Eq = K(Q) + U c ff |(T.;)| 2 + 
U e g/4:. As shown in Fig. |3 an excellent mapping results 
from t = -0.6 eV, U cS = 1.72 eV, A JT = gQ* = 0.87 eV, 
and E z = 0.12 eV. Remarkably, even 9 T = tan" 1 T x /T z 
obtained from the effective Hamiltonian also agrees well 
with the LDA+C/ results. That is, the effective Hamilto- 
nian accurately reproduce not only all the LDA+C7 en- 
ergies, but also the wavefunctions. 

Further insight into the relative importance of e-e and 
e-1 interactions now can be obtained by examining the 
contribution of each term of the effective Hamiltonian, 
Eq. (J2J) , to the HF total energy, as shown in Table [I] 
First, consider Ao (first row), the energy gain purely due 
to the OO formation (finite (Tj)) in the absence of JT 
distortion. Consistent with the above LDA+C/ results, 
e-e interaction alone (—173 meV) overwhelms the kinetic 
energy cost (132 meV) and induces OO with the total en- 
ergy gain of —50 meV. More intriguingly, the additional 
energy gain via cooperative JT distortion Q*, referred to 
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TABLE I: Contribution of the energy terms in Eq. J5J to A 
and Aq* (see text), in unit of meV per formula unit. 





Total 


u eS 
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E x 


t 




-50 


-173 








-9 


132 


A Q . 


-215 


-129 


-356 


96 


16 


158 



as Aq* (second row), consists of considerable further con- 
tribution from e-e interaction (—129 meV). Without the 
e-e interaction, the JT coupling (—356 meV) alone is in- 
sufficient to overcome the energy cost of lattice (96 meV) 
and kinetic (290 meV) energy. That is, c-c interaction, 
whose contribution is enhanced due to increased (Tj), 
serves as a hidden driving force for octahedral distortion. 
Overall, U e s is more important to 00 in LaMn03 than 
the JT coupling, while their effects on cooperative octa- 
hedral distortion are comparable. 

The present results place stringent constraints on any 
realistic theories of the manganites and interpretations of 
their excitation spectra, for example, current controversy 
on the recently observed 120 — 160 meV Raman shifts of 
incident photons resonant at E rcs ~ 2 eV 
the one hand, the excitations were interpreted as orbital 
waves, or 'orbitons', derived from the pseudo-spin su- 
perexchange model in the large E/ G ff limit of Eq. (El |2fJ . 
With the orbiton spectrum gap ~ 2.5J or b + Ajt U Lii| 
and superexchange coupling constant J or b — 40 — 50 meV 
(oc t 2 / Ue.f r ), this scenario requires small Ajt < 50 meV 
|2l La, |30( . On the other hand, in the JT scenario the 
Raman shifts were attributed to two-phonon processes 
(single phonon frequency ~ 60 — 80 meV) induced by 
a phonon assisted on-sit e d£ — ► df transition, requiring 
large Ajt — E rcs ~ 2 eV [l7|,|31|. In contrast, our quan- 
titative results (Ajt — 0.9 eV and U c g — 1.7 eV) show 
that the spin-majority e g states are in the intermediate e- 
e interaction regime with comparable JT coupling. Thus, 
a more reasonable picture would be two-phonon processes 
mediated by optically active inter-site djdlj — ► df dfc tran- 
sition with E xes ~ 2 eV |l4l Il5j . Direct experimental ver- 
ification of our results includes probing, e.g., on-site d-d 
transition (~ 0.9 eV) via inelastic X-ray scattering, or 8? 
via nuclear magnetic resonance. 

In summary, we have quantified the relative impor- 
tance of e-e and e-1 interactions in ordering orbitals in 
LaMnOa using a new theoretical approach. A realis- 
tic effective Hamiltonian resulting from this quantita- 
tive method reproduces consistently both LDA+L 7 en- 
ergies and wavefunctions. Intermediate c-e interaction 
{UeS — 1-7 eV) is found to play a crucial role in inducing 
OO and localizing the electrons, which in turn enhances 
the JT interaction (Ajt — 0.9 eV) and stabilizes JT dis- 
tortion. Furthermore, a clear "signature" of competition 
between e-e and e-1 interactions is given via the orbital 
mixing angle 8t < Oq. Experimental means to directly 
clarify the relative strengths of the leading mechanisms 



are suggested. The developed general theoretical scheme 
can be applied to other strongly correlated materials. 

We are grateful to P. Allen, E. Dagotto, S. Grenier, J. 
Hill, A.J. Millis, A. Moreo, G. Sawatzky, D. Singh, and 
J. Thomas for helpful discussions. W.Y. thanks M. Lu- 
faso and P. Woodward for providing SPuDS [2(| and T. 
Chatterji for providing structural data [13] ■ Brookhaven 
National Laboratory is supported by U.S. Department of 
Energy under Contract No. DE-AC02-98CH1-886. This 
work is partially supported by DOE-CMSN. 



[1] E. Dagotto, Nanoscale Phase Separation and Colossal 
Magnetoresistance (Springer Series in Solid State Sci- 
ences vol. 136, Springer, 2003). 

[2] T. Hotta, A. L. Malvezzi, and E. Dagotto, Phys. Rev. B 
62, 9432 (2000). 

[3] W.-G. Yin, H. Q. Lin, and C. D. Gong, Phys. Rev. Lett. 
87, 047204 (2001). 

[4] J. Bala and A. M. Oles, Phys. Rev. B 62, R6085 (2000). 

[5] M. V. Mostovoy and D. I. Khomskii, Phys. Rev. Lett. 
92, 167201 (2004). 

[6] K. H. Ahn and A. J. Millis, Phys. Rev. B 61, 13545 
(2000). 

[7] S. Okamoto, S. Ishihara, and S. Maekawa, Phys. Rev. B 

65, 144403 (2002). 
[8] R. Tyer, W. M. Temmerman, Z. Szotek, G. Banach, 

A. Svane, L. Petit, and G. A. Gehring, Europhys. Lett. 
65, 519 (2004). 

[9] A. J. Millis, Nature (London) 392, 147 (1998). 
[10] C. M. Varma, Phys. Rev. B 54, 7328 (1996). 
[11] A. S. Alexandrov, and A. M. Bratkovsky, Phys. Rev. 

Lett. 84, 2043 (2000). 
[12] T. V. Ramakrishnan, H. R. Krishnamurthy, S. R. Hassan, 

and G. Venketeswara Pai, Phys. Rev. Lett. 92, 157203 

(2004). 

[13] Y. Murakami, J. P. Hill, D. Gibbs, M. Blume, I. Koyama, 
M. Tanaka, H. Kawata, T. Arima, Y. Tokura, K. Hirota, 
et al., Phys. Rev. Lett. 81, 582 (1998). 

[14] S. Grenier, J. P. Hill, V. Kiryukhin, W. Ku, Y.-J. Kim, 
K. J. Thomas, S.-W. Cheong, Y. Tokura, Y. Tomioka, 
D. Casa, et al., Phys. Rev. Lett. 94, 047203 (2005). 

[15] N. N. Kovaleva, A. V. Boris, C. Bernhard, A. Ku- 
lakov, A. Pimenov, A. M. Balbashov, G. Khaliullin, and 

B. Keimer, Phys. Rev. Lett. 93, 147204 (2004). 

[16] R. Rauer, M. Riibhausen, and K. Dorr, to be published. 

[17] R. Kriiger, B. Schulz, S. Naler, R. Rauer, D. Budelmann, 
J. Backstrom, K. H. Kim, S.-W. Cheong, V. Perebeinos, 
and M. Riibhausen, Phys. Rev. Lett. 92, 097203 (2004). 

[18] E. Saitoh, S. Okamoto, K. T. Takahashi, K. Tobe, 
K. Yamamoto, T. Kimura, S. Ishihara, S. Maekawa, and 
Y. Tokura, Nature (London) 410, 180 (2001). 

[19] J. Kanamori, J. Appl. Phys. 31, 14S (1960). 

[20] K. Kugel and D. Khomskii, Sov. Phys. JETP 37, 725 
(1973). 

[21] W. Ku, H. Rosner, W. E. Pickett, and R. T. Scalettar, 

Phys. Rev. Lett. 89, 167204 (2002). 
[22] V. I. Anisimov, F. Aryassetiawan, and A. I. Lichtenstein, 

J. Phys.: Condens. Matter 9, 767 (1997). 
[23] LDA+f/ (U=8 eV, J=0.11U) was shown to appropriately 



5 



describe the (ir, tt, 0) orbital and (0, 0, n) magnetic struc- 
tures of LaMnOs, see I. S. Elfimov and V. I. Anisimov 
and G. A. Sawatzky, Phys. Rev. Lett. 82, 4264 (1999). 

[24] All our LDA+U results were obtained with the linearized 
augmented plane wave method I25I for the A-type anti- 
ferromagnetic phase with experimental volume and fixed 
J=0.11U. Hypothetic crystal structures were prepared 
with SPuDS IS. 

[25] P. Blaha et ai, Comput. Phys. Commun. 147, 71 (2002). 

[26] M. W. Lufaso and P. M. Woodward, Acta Cryst. B 60, 



10 (2004). 

[27] T. Chatterji, F. Fauth, B. Ouladdiaf, P. Mandal, and 
B. Ghosh, Phys. Rev. B 68, 052406 (2003). 

[28] W. E. Pickett and D. J. Singh, Phys. Rev. B 53, 1146 
(1996). 

[29] W. Ku et al., to be published. 

[30] J. van den Brink, Phys. Rev. Lett. 87, 217202 (2001). 
[31] P. B. Allen and V. Perebeinos, Phys. Rev. Lett. 83, 4828 
(1999). 



